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ABSTRACT 

Bayesian model averaging is a procedure to obtain parameter constraints that account 
for the uncertainty about the correct cosmological model. We use recent cosmological 
observations and Bayesian model averaging to derive tight limits on the curvature 
parameter, as well as robust lower bounds on the curvature radius of the Universe 
and its minimum size, while allowing for the possibility of an evolving dark energy 
component. Because flat models are favoured by Bayesian model selection, we find 
that model-averaged constraints on the curvature and size of the Universe can be 
considerably stronger than non model-averaged ones. For the most conservative prior 
choice (based on inflationary considerations), our procedure improves on non model- 
averaged constraints on the curvature by a factor of ~ 2. The curvature scale of the 
Universe is conservatively constrained to be R c > 42 Gpc (99%), corresponding to a 
lower limit to the number of Hubble spheres in the Universe Njj > 251 (99%). 
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1 INTRODUCTION 

One of the most spectacular advancements of observational 
cosmology over the past two decades has been the ability 
to measure the spatial geometry of the universe with un- 
precedented accuracy. In the Friedmann-Robertson- Walker 
(FRW) Universe there are only three discrete possibilities for 
the underlying geometry, namely flat, open or closed. The 
amount of curvature is usually characterized by the curva- 
ture parameter Sl K : if S} K < the geometry of spatial sections 
is spherical (i.e., the Universe is closed) and the Universe has 
a finite size. If instead il K > the geometry is hyperbolic 
(i.e., the Universe is open), while for £l K = spatial sections 
are flat. In both the two latter cases, the spatial extent of 
the Universe is infinite. 

Limits on the value of Q, K can be derived in a geo- 
metrical way by observing the angular size subtended by 
cosmological features of known physical length, such as 
the acoustic peaks in the cosmic microwave background 
(CMB) and the corresponding baryonic acoustic oscilla- 
tions (BAO) in the distribution of large scale structures. 
Furthermore, type la supernovae (SNIa) can be used as 
standard candles to determine the luminosity distance as a 
function of redshift. A combination of these three probes 
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has been succesfully used to set very tight limits to the 
curvature parameter, which is now constrained at better 
than the ~ 10 -3 level. For example, iKomatsu et al.l (|2011f ) 
find tt K = -0.0057lo.oo68 at 68 % CL, employ ing a com- 
binat ion of WMAP7, B AO (|Reid et all l20ld ) and SNIa 
data jHicken et alj|2009h . 

Impressive as such limits are, they assume the Universe 
to be curved, and carry out parameter inference on the quan- 
tity describing curvature. A different methodological per- 
spective is required to go beyond that assumption: Bayesian 
model comparison can be used to obt ain the posterior 
prob ability of the Universe being curved jVardanvan et all 
120091 ). The purpose o f this paper is to furthe r expand the 
approach adopted in IVardanvan et all (|2009l ) by deriving 
model-averaged limits on the curvature of the Universe, 
which fully account for the uncertainty in selecting the cor- 
rect model for the FRW Universe. Given current data, flat 
models are preferred by Bayesian model selection from an 
Occam's razor perspective, and therefore most of the proba- 
bility mass becomes concentrated in models with vanishing 
spatial curvature. As we demonstrate, model-averaged lim- 
its on the curvature parameter are tighter than non model- 
averaged constraints, because of this "concentration of prob- 
ability" effect. 

This paper is organized as follows: in section [5] we out- 
line the Bayesian model averaging framework. We describe 
the data sets used for the analysis and the models included 
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in section [3] We then present our results in section [4j where 
we also discuss the dependency on prior choice. We give our 
conclusions in section [5] 



is given by 



P (6\d) = ^2p(9\d,Mi) P (Mi\d) 



(5) 



P (Mi\d)Y,Bii P (8\d,Mi), 



(6) 



2 BAYESIAN MODEL AVERAGING 

If one is interested in determining the parameters 6 of a 
given model Ai, then the relevant quantity is the posterior 
distribution for 6 under that model, which is given by Bayes 
theorem as 



p(6\d,M) = 



p(d\0,M)p(6\M) 
P (d\M) 



(1) 



where the explicit conditioning on Ai indicates that the pos- 
terior pdf for 9 given data d, p(6\d, Ai), is conditional on 
having assumed a specific model Ai. This is the usual pa- 
rameter inference step, and often the first level of inference in 
a problem (i.e., determining the constraints on the model's 
parameters). 

The second level of inference is Bayesian model com- 
parison, which aims to determine the relative probability of 
models themselves. The posterior probability of a model Aii 
given the data, p(Aii\d) is related to the Bayesian evidence 
(or model likelihood) p(d\Aii) by 



p(Mi\d) = 



p(d\Mi)p(Mi) 
p(d) 



(2) 



where p(Aii) is the prior for model Aii, p(d) = 
£V p(d\ Mi)p(Mi) is a normalization constant (where the 
sum runs over all available models) and 



p(d\Mi)= dep{d\e,M t )p{8\M t ) 



(3) 



is the Bayesian evidence, which appears as a normalization 
factor in Eq. JJJ. Given two competing models Aii, Aij, the 
change in their relative probability in going from the prior 
(i.e, before we see the data) to the model posterior (after 
the data have been taken into account via the likelihood) is 
given by the Bayes factor Bij : 



B i 



p(d\Mi) 
p(d\M s ) 



(4) 



where large (small) values of Bij denote a preference for 
Aii (Aij). The 'Jeffreys' scale' gives an empirical scale 
for translating the values of ln_B;j into strengths of belief, 
with thresholds |lnBy| = 1.0,2.5,5.0 separating levels of 
inconclusive, w eak, moderate and strong evidence, respec- 
tively, see e.g. iTrottal |2008l ). Recently, the framework of 
model comparison has been extende d to include the possi- 
bility of 'unknown models' discovery (|Starkman et alj|2008l ; 
iMarch et afll201Cl ). 

The third level of inference is represented by Bayesian 
model averaging (BMA) , whose purpose is to determine con- 
straints on common parameters among the models being 
considered (Aii, with i = 1, . . . , N) accounting for the un- 
certainty in selecting the correct model. This is the most 
general inference one can obtain on the parameters values 
(at least as long as the list of models is reasonably complete). 
The model-averaged posterior distribution for parameter 6 



where in the second equality we have replaced the models 
posterior probabilities, p(Aii\d) by the Bayes factors with 
respect to a reference model, here Aii, and further assumed 
that the prior probabilities for the iV models are all equal, 
i.e. p(Aii) — 1/N, (i — 1,...,N). With this assumption, 
the posterior for Aii is given by 

p(Mi\d)= * . (7) 

BMA h as been appl i ed to the dark energy equation of 
st ate in iLiddle et al.l (I2006T). to the scalar spectral index 
in IParkinson fc Liddlel ( 2010h and in the c ontext of weak 
lensing and Sunyaev-Zel'dovich effect data in Marshall et al. 
(|2003h . 

2.1 Computation of the Bayes factors 

Given two or more models, computing the Bayes fac- 
tors entering Eq. (0 requires the evaluation of the multi- 
dimensional integral of Eq. ([3]). Several algorithms are avail- 
able today to compute the Bayesian evidence numerically. 
Here we are interested in the case where the models are 
nested within each other, i.e. where one of the models is ob- 
tained from a more complicated one for a speficic choice of 
some of the parameters of the latter. In our case, the extra 
parameters are the curvature, £l K , and/or the dark energy 
equation of state parameters, wo or w a , depending on the 
model under consideration. For example, a curved Universe 
reverts to a flat one for fi K = 0, or an evolving dark energy 
equation of state reverts to a cosmological constant model 
for wo = —l,w a — 0. In this case, the Bayes factor between 
models Aii and Aij can be written in all generality as 

p(-3\d,Mj) 



B, 



p{0\Mj 



(8) 



where we have split the more complicated model's parame- 
ters as 6 — (-0, i9), and i? are the extra parameters of model 
Aij, which reduces to the simpler model Aii for # = 
This express ion is known as the Savage-Dick ey density ratio 
(SDDR, see IVerdinelli fc Wassermarl (| 19951 ) and references 
therein. For cosmological applications, see ITrottal |2007h ). 
The numerator is simply the marginal posterior for i9, evalu- 
ated at the value, $ = (which can easily be obtained with 
standard Markov Chain Monte Carlo techniques), while the 
denominator is the prior density for the extra parameters $ 
under the more complicated model, evaluated at the same 
point. 

Once the Bayes factor for nested models which differ by 
one parameter at the time has been obtained using Eq. ([S]), 
the Bayes factor between other models which have two or 
more nested parameters between them can be easily derived. 
If model Aii has one more parameter than model Aik, which 
in turn has one more parameter than Aij , the Bayes factor 
between models i and j is given by 

/<•,, /;./ X Ik .- (9) 
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where the Bayes factor B^ and B^j can be obta ined via the 
SDDR . A similar technique has been adopted in lKunz et al.l 
l|200eh . 



3 MODELS AND DATA SETS 
3.1 Models and priors 

We work in the framework of FRW cosmologies including 
a cold dark matter (CDM) component, a possible curva- 
ture and a dark energy component. We will consider open 
(k = —1) and closed (k = 1) Universes as two separate 
models, as this allows us to adopt a prior on the curvature 
parameter Q K which is uniform in log(Sl K ), as explained be- 
low. A model is fully specified by the choice of curvature and 
dark energy parameterization, alongside their respective pri- 
ors. We consider the possibilities listed in Table [1] a model 
is defined by picking a choice from the upper part of the ta- 
ble (curvature) and one from the lower part (dark energy). 
From now on, a single model will be denoted by a pair of 
labels in the subscript, each referring to the prior choice for 
the curvature and dark energy sector. So, for example, the 
model Mfa denotes a flat Universe (F) with a cosmolog- 
ical constant (A). Correspondingly, Bayes factors between 
two models will have two pairs of labels in the subscript, 
separated by a comma for clarity, e.g. Bow, fa denotes the 
Bayes factor between Mow and Mfa- 

For non-flat models, we consider two different prior 
choices for the curvature parameter: a uniform prior in the 
range —1 ^ fi K ^ 1 (the 'Astronomer's prior') and a uniform 
prior in the range —5 ^ log |fi K | ^ (the 'Curvature scale 
prior'). The Astronomer's prior is motivated by basic con- 
sistency with observable properties of the Universe, such as 
the age of the oldest objects, while the Cur vature scale prior 
is bas ed on an inflationary scenario, see IVardanvan et all 
for full details. For the dark energy equation of state, 
we adopt the parameterization 



Csee lScott fc Zibinl < |2006f > for details), 



w(z) = Wo + w a z/(l + z), 



(10) 



with two free parameters, wo,w a . Model W has w a = and 
a uniform prior on rao in the range given in Table [TJ while 
model W Z further allows for w a 7^ 0. The prior ranges for 
wo, w a are sufficiently wide to enclose the support of the like- 
lihood function, but not too large in order to avoid a very 
strong Occam's razor effect against evolving dark energy 
models. Finally, parameters that are common to all mod- 
els (such as the amplitude of primordial fluctuations or the 
baryonic density) are irrelevant for the model comparison, 
as shown by Eq. JSJ, and therefore the choice of priors on 
them is unproblematic. For each choice of curvature prior, 
we thus consider a total of 9 cosmological models. 

In the following, we will derive model-averaged con- 
straints on the curvature parameter, Q. K , the curvature ra- 
dius R c , given by 



Rc = 



Ho |n K |V2 



(11) 



(where c is the speed of light and Ho the Hubble constant 
today in km/s/Mpc) and the number of Hubble spheres in 
the Universe, Nu, defined as the ratio of the present vol- 
ume of the spatial slice to the apparent particle horizon 



Nn 



2tt 



2 X - sin(2 X ) ' 



(12) 



where \ is the comoving radial distance. 

In some models, the value of some of the parameters is 
fixed. E.g., for flat models the curvature parameter vanishes, 
£l K — (for the Astronomer's prior), while for open and flat 
model R c , Nu — > 00. In such cases, the posterior probability 
mass associated with that model in the model-averaged ex- 
pression gets concetrated in a Dirac delta function 5 around 
the fixed value of the parameter. So in the case of curvature, 
for example, the BMA expression (JSJ) becomes: 



p(n„|d) = p (Mfa\cI)\ 1 Yl p(n J ,|d,M y )B ij,FA + 

[i^F,j=A,W,WZ] 

Fj,FA\, 

[j=A,W,WZ] 

(13) 

where the fixed value £1* = for the Astronomer's prior 
and Q% = 10 _J for the curvature scale prior. In the above 
expression we have taken the flat, ACDM model FA as our 
reference model, and computed all Bayes factors with re- 
spect to it. 



3.2 Data sets employed 

The angular position of the first acoustic peak in the CMB 
power spectrum and the same acoustic signature in the cor- 
relation function of the galaxies provide us with standard 
rulers at high and low redshifts, respectively. The measure- 
ments of the standard rulers in the direction perpendicular 
to the line of sight are used to constrain the angular di- 
ameter distance Da{z) = (1 + z)~ 2 Dl(z), where Dl is the 
luminosity distance, given by 



Dl(z) = 



= (1 + 2) sin 



dx 
H(x) 



(14) 



The Hubble function H(z) is expressed in terms of the 
present-day matter-energy content of the Universe as fol- 
lows: 



H 2 (z) = (fi m (l + zf + n r (l + zf + Q K (1 + zf 



+ Ode exp ( 3 



z 1 + w(x) 
1 + x 



dx 



(15) 



The dark energy time evolution is described by the present- 
day dark energy density in units of the critical density, Qde, 
and by its equation of state, w{z), as given in Eq. (|10|l . In 
extracting constraints on cosmological parameters from lu- 
minosity or angular diameter distance measurements, one 
has to be careful to consider the potential impact of de- 
generacies between the assumed models. In this case, the 
strong degener acy between cur v ature and dark energy evo- 
lution (see e.g. IClarkson et al.l f2007])) is at least partially 
accounted for by admitting in our space of models an evolv- 
ing dark energy equation of state. 

W e include the WMAP 5-year data (|Dunklev et al.l 
I2009D via their constraints on the shift parameters 
an d the baryon density , following the method employed 
in lKomatsu et al.l (|2009n . We notice that adopting WMAP 
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Model 


Parameters and priors 




Curvature sector 
Astronomer's prior Curvature scale prior 


F (flat) 
O (open) 
C (closed) 


n K = o n K = icr b 

^ Q K ^ 1 (uniform) —5 logf2 K ^ (uniform) 
-1 < n K ^ (uniform) -5 sj -logQ^ < (uniform) 




Dark energy sector 


A 

W 

WZ 


u>o = —l,w a = 

—2 ^ u>o ^ —1/3 (uniform), w a = 

—2 ^ uiq ^ —1/3 (uniform), —1.33 ^ w a 1.33 (uniform) 



Table 1. Prior choices for the curvature (top half) and dark energy (bottom half) parameters considered in the analysis. A model is 
fully specified by selecting a prior choice from the top and one from the bottom of the table, thus defining both the curvature and the 
dark energy sectors. 



7-year data is not expected to change significantly our re- 
sults, because of the fundamental geometrical degeneracy in 
the CMB. We make use of the SDSS and 2dFGRS baryoni c 
acoustic scale measurements following lPercival et al.l (|2007i ). 
The scale of BAOs is used to estimate the quantity 

Dy(z) = (cz bao (l + zf^^f) , (16) 

for Zbao = 0.2 and Zb ao = 0.35. We use two Gaus- 
sian data points with mean value and standard devia- 
tion r s /Dv(0.2) = 0.1980 ± 0.005 8 and r s /D v (0.35) = 
0.1094 ±0.0033 (|Percival et alj|2007j ). where r 3 is the acous- 
tic sound horizon. We also add the Hubble Key Project 
determination of the Hubble constant today, as a Gaus- 
sian datum with mean and standard deviation Ho = 
72 ± 8 km/s/Mpc (|Freedman et al.ll200ll ). SN type la data 
are included in the form of the UNION08 data set sam- 
ple jKowalski et alj|2008h . 

We employ a Metropolis-Hastings Markov Chain Monte 
Carlo procedure to derive the posterior distribution for the 
parameters in our model, and to compute the Bayes factors 
necessary for BMA via Eq. ([H|. We take flat priors on the 
following quantities: fi m /i , Qbh 2 , wo, w a (whenever wo, w a 
is not fixed to —1, 0, respectively). The prior bounds on the 
first 2 parameters are irrelevant, as the posterior is well con- 
strained within the prior. The Bayes factors are obtained as 
the mean of the Bayes factors obtained from of 8 indepen- 
dent reconstructions with MCMC, while their uncertainty is 
estimated from the variance of the values from the 8 runs. 



4 RESULTS AND DISCUSSION 

Table [5] gives the results of our model comparison between 
the flat ACDM model and the 8 alternative models con- 
sidered, for both choices of curvature prior. As all values 
of ln_B > 0, the flat ACDM model is the one preferred 
by the data, as expected. The strength of evidence in its 
favour depends on the chosen prior, with the curvature scale 
prior giving in general a smaller evidence against the alter- 
native models, since its Occam's razor effect on curvature 
is smalle r than for the As t ronom er's prior, as discussed in 
detail in IVardanvan et~aL | l|2009D . For example, when com- 
paring a flat ACDM model with a closed ACDM Universe, 
the odds in favour of the former are ~ 100 : 1 under the 
Astronomer's prior (moderate evidence), while only 2 : 1 
under the more conservative Curvature scale prior (inconclu- 
sive evidence). We also observe weak to moderate evidence 



against the evolving dark energy models as compared to 
models with a cosmological constant (both flat and curved). 
The posterior probability for fiat models is given by 

P(n = 0\d) = ( 1+ W W ^a V' 

which gives p(n = 0|d) = 0.986 for the Astronomer's prior, 
starting from a prior probability p(n = 0) = 3/9 « 0.33 
(as there are 3 flat models among the 9 models we con- 
sider here). For the more conservative Curvature scale prior, 
we obtain a posterior probability for flatness of only p(n = 
0\d) = 0.462. The origin of this difference can once more 
be traced back to the less pronounced Occam's razor effect 
against non-flat models for the Curvature scale prior. 

Using the values of the Bayes factors given in Table [21 
we derive in Table [3] model-averaged constraints on the 
quantities of interest (Q K ,Nu,Rc)- The table also gives the 
non-model averaged constraints for the most conservative 
model (with a possibly time-evolving dark energy equation 
of state) for comparison. The model-averaged posterior for 
the curvature parameter is also plotted in Fig. Q] for both 
choices of priors, where the spike in both panels represents 
the posterior probability mass concentrated in the flat mod- 
els (~ 99% for the Astronomer's prior and ~ 46% for the 
Curvature scale prior). Because the posterior probability of 
fiat models gets concentrated in a delta-function at £7 K = 
(for the Astronomer's prior) or |n K | = 10 -5 (for the Curva- 
ture scale prior), the model-averaged constraints on f2« can 
produce tighter intervals than the usual results obtained un- 
der the assumption that fi K 7^ 0. This "concentration of 
probability" onto the simpler model is a consequence of the 
Occam's razor effect implicit in Bayesian model selection. 
A similar effect has been observed in model-a veraged con- 
strain ts on the dark energy equation of state (|Liddle et al.l 
120061 ). Because more than 95% of the posterior probabil- 
ity is concentrated in flat models under the Astronomer's 
prior, 95% limits on the curvature parameter are not de- 
fined for this prior. At 99%, we find that \Q K \ < 2 x 10~ 4 , 
while the number of Hubble spheres is greater than 398 and 
the radius of curvature larger than 68 Gpc (all at 99%). The 
Curvarture scale prior is more conservative, in that it penal- 
izes less strongly non-flat models. As a consequence, we find 
that -0.9 x 10~ 2 < D. K < 1.0 x 1.0~ 2 (99% model-averaged 
region), while the number of Hubble spheres is Nu > 251 
and the curvature radius R c > 42 Gpc. Even under this 
more conservative prior, the model-averaged constraint on 
curvature is a factor ~ 2 better than the usual, non model- 
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Prior 


CWZ 


CW 


CA 


OWZ 


OW 


OA 


FWZ 


FW 


p(K = 0\d) 


Astronomer's 
Curvature 


8.28 ± 0.09 
3.28 ± 0.1 


6.60 ±0.09 
2.79 ±0.1 


4.61 ±0.09 
0.73 ±0.04 


7.41 ±0.06 
2.82 ±0.05 


7.07 ±0.05 
2.31 ±0.04 


5.63 ± 0.03 
0.44 ± 0.04 


2.67 ±0.03 
2.67 ±0.03 


2.13 ±0.03 
2.13 ±0.03 


0.986 ± 0.005 
0.462 ± 0.006 



Table 2. Log of the Bayes factors (lnB) between the flat ACDM model and alternative models with curvature and/or an evolving dark 
energy equation of state (indicated by the label in the top line), for both our choices of priors for the curvature sector. The last column 
gives the posterior probability of the Universe being flat. 




-0.02 -0.01 0.01 io- 6 io-« io- 3 



Figure 1. Model-averaged posterior probability distribution for 
the curvature parameter, including all 9 models considered in the 
analysis, assuming the Astronomers' prior (left panel) and the 
Curvature scale prior (right panel) for the curvature parameter. In 
the right panel, the solid line applies to closed Universes (Q K < 0), 
while the dotted line to open Universes (fi K > 0). The peaks 
represent the Dirac delta function encompassing the probability 
mass associated with flat models. 



Quantity 


Astronomers' prior 
95% 99% 


Curvature scale prior 
95% 99% 


Model-averaged constraints 


n K X IO" 3 
Rc (Gpc) 
N v 


N/A [-0.2,0.2] 
N/A > 68 
N/A > 398 


[-4.4, 2.0] [-8.9, 9.9] 

> 88 > 42 

> 1000 > 251 


Non model-averaged constraints (WZ model) 


n K X io -3 

Rc (Gpc) 
Nu 


[-14,9] [-19,17] 

> 34 > 32 

> 33 > 29 


[-9, 7] [-14, 13] 

> 30 > 28 

> 38 > 32 



Table 3. Top section: model-averaged parameter constraints for 
the curvature parameter, the curvature scale and the number of 
Hubble spheres in the Universe. For the Astronomer's prior, ~ 
97% of the posterior probability is concentrated in the flat models, 
and therefore the 95% limit for the model-averaged parameters 
is not defined. Bottom section: non model-averaged constraints 
for the most conservative case (evolving dark energy model), for 
comparison. 

averaged results, which (from a similar colle ction of data) is 
appro ximately \£l K \ ^ 1.7 X 10~ 2 at 99% (|Komatsu et al.l 
|2011| ). The model -averaged constraint on the number of 
Hubble spheres is a factor ~ 40 stronger than the non model- 
averaged on e: in the latter c a se, an d using the Curvature 
scale prior, IVardanvan et al.1 l|2009l ) found N v > 6.2 (at 
99%). 



5 CONCLUSIONS 

We have applied the formalism of Bayesian model averag- 
ing to the problem of constraining the curvature and size of 
the Universe. By employing the Savage-Dickey density ra- 



tio, we have obtained model-averaged constraints at almost 
no additional computational effort than what is needed for 
parameter estimation. We have demonstrated how model- 
averaged constraints on the curvature and minimum size of 
the Universe can be considerably tigther than non model- 
averaged ones. This is a consequence of the fact that flat 
models are preferred by Bayesian model selection, although 
the strength of such preference is fairly strongly dependent 
on the choice of prior for the curvature parameter. 

We have considered two classes of priors that are 
based on physical and theoretical considerations. We found 
that even the most conservative prior choice gives model- 
averaged constraints on curvature that are a factor of ~ 2 
better than non model-averaged intervals. A more aggressive 
prior choice (the Astronomer's prior) leads to an improve- 
ment in the constraints on Q K by a factor ~ 100, giving 
|fi„| < 2 x 10~ 4 at 99%. The minimum size of the Uni- 
verse is robustly constrained to encompass Nu > 251 Hub- 
ble spheres, an improvement of a factor ~ 40 on previous 
constraints. Finally, the radius of curvature of spatial section 
is found to be R c > 42 Gpc. 

Bayesian model averaging gives the most general pa- 
rameter constraints, which fully account for the uncertainty 
in the selection of the correct underlying cosmological model. 
It remains imperative (like in any good Bayesian analysis) 
to study the dependency of the results of the chosen priors, 
which are more important in Bayesian model selection (and 
model averaging) than they are in the usual parameter in- 
ference framework. We believe that the formalism presented 
here can be employed successfully in a large variety of cos- 
mological problems. 
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